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ABSTRACT 

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred ob- 
servational properties, rely on the assumption that angular momentum transport is opposite to the radial angular 
frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by 
construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrody- 
namic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational 
instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency in- 
creases with radius. In order to shed light into physically viable mechanisms for angular momentum transport 
in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differen- 
tially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We 
isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evo- 
lution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be 
amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism 
to transport significant angular momentum (inward). These findings are consistent with the results obtained 
in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of 
efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of 
turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard 
framework of turbulent shear viscosity. 

Subject headings: accretion, accretion disks — instabilities — magnetohydrodynamics (MHD) — turbulence 



1. INTRODUCTION 

Basic arguments suggest that the angular frequency Q.(r) 
of an accretion disk surrounding a weakly magnetized star 
must attain a maximum value, Qmax = ^(Xh ), and de- 
creas e inward (or at least remain constant; see IMedvedevI 
120041) to match the angular f requency of t he star at the stel- 
lar rad i us flvf(r»); see, e .g., IFrank et al] ( 120021) : IHartmannl 
(120091) : lArmitagj (120101) . The inner disk region, where 
r < rb and dQ./dr > 0, is referred to as the accre- 
tion disk boun dary layer. Standard accretion disk theory 
([Shakura & Su nvaevI 119731) predicts that half of the energy 
released in the accretion process takes place in this region, 
estimated to be a fraction of the stellar radius. The spec- 
trum of the radi ated energy depends on the detailed properties 
of this layer (Naravan & PoDham"1993"; 'Popham & NarayanI 
■ [T995; Popham etai. 1996: Popham & Sunyaev 2001): thus 
understanding the various processes that determine its prop- 
erties (see, e.g., Piro & Bildsten 2004: IBalsara et all 120091: 
llnogamov & Sunyaevi201 0) is of fundamental importance. 

Most detailed calculations for determining the structure 
of the boundary layer rely on effective models for turbu- 
lent angular momentum transport. These models are usu- 
ally built as a turbulent version of the Newtonian viscous 
stress bet ween fluid layer s in a differentially rotating lami- 
nar flow (ILandau & Lifsh itz 1959), and thus assume a lin- 
ear relationship between the stress and the angular fre- 
quency gradient (Lynden-Bell & Pringle 1974). This assump- 
tion, however, seems at odds with the properties of mag- 
netohydrodynamic (MHD) turbulence revealed by numeri- 
cal simulat ions of accretion d isks ( Abramowicz et al. 199^ 
lArmitage 2002t ISteinacker & Papaloizou 2002: Pessah et al] 
12008) which show that angular momentum transport is inef- 
ficient in regions of the disk where dQ./dr > 0, which are 
stable to the standard magnetorotational instability (MRI: see 



IBalbus & HawlevlfT99TlfT998h . 

Motivated by the need of a deeper understanding of the be- 
havior of an MHD fluid in a differentially rotating background 
that deviates from a Keplerian profile, we study the dynamics 
of MHD waves in configurations that are stable to the standard 
MRI. Employing the shearing-sheet framework, we show that 
transient amplification of shearing MHD waves can generate 
magnetic energy without leading to a substantial generation of 
hydromagnetic stresses. We discuss the implications of these 
findings. 

2. ASSUMPTIONS AND LOCAL MODEL FOR MHD DISK 

We focus our attention on a subsonic, weakly magne- 
tized fluid for which the ram pressure and magnetic pres- 
sure remain small compared to their thermal counterpart. As 
a first approximation, we thus consider a differentially ro- 
tating fluid with angular frequency £1 - Q(r)z and con- 
stant background density po. This is a reasonable as sump- 
tion in light of the results presented by Armitage (2002) and 
Steinacker & Papaloizou (2002), who carried out numerical 
simulations of boundary layers of unstratified accretion disks 
and found that the density fluctuations throughout the simula- 
tions are in general quite small. 

We work in the framework of the shearing-sheet 
approximation dffilil 11878 : Goldreic h & Tremaind 119781: 
INarav an et al.lll98"7l : IHawlev et aljl995i) . where the equations 
describing an incompressible MHD fluid in a corotating frame 
are given by 

9 VP 

d,v + (v ■ V)v - -2QoZ X V + 2qQ.^yXX 



47rpo 

d,B + (v ■ V)B = (B ■ V)v + rjV^B . (2) 
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Here, v{x, t) and B{x, t), with V • v = V ■ B = 0, stand for 
the velocity and magnetic fields; v and rj denote the kinematic 
viscosity and resistivity; and Qq = Q(ro) is the corotating an- 
gular frequency at a fiducial radius kq. The first and second 
terms on the right hand side of equation ([T]) correspond to the 
Coriolis and tidal forces, respectively. The local pressure P 
can be found by the divergence-less condition of the velocity 
field. Recalling that the local density po is assumed to be con- 
stant, and in order to simplify notation, hereafter we redefine 
the symbols denoting the pressure and magnetic field in such 
a way that and P/po -> P and B/(47rpo)'^^ B. 
We decompose the flow into mean and fluctuations as 



v{x,t) = Ui{x) + u{x,f), 
B(x,t) = Bo(t) + b(x,t). 

The leading order background velocity is Ui{x) = 
where the shear parameter q is given by 



£/lnQ 



q = - 



dlnr 



(3) 
(4) 

-q Qoxy, 



(5) 



The homogeneous background magnetic field is in general 
a function of time and it evolves according to the induction 
equation i.e., d/Bo - -^QoBo.^j. 

The substitution of equations Q and (HI into ([T]i and (|2|i 
leads to a non-linear system for the dynamical evolution of 
the perturbations u(x, t) and b{x, t). However, as pointed out 
in Goodma n & Xui (|1994), all the non-linear terms in the re- 
sulting equations vanish identically if we consider the evolu- 
tion of a single Fourier modf]. In this case we obtain 

D,u - q Q.ouJ ^ Bo-Vb + vV^H - 2QoZ xu-VP, (6) 

D,b + q Q.obJ = Bo ■ Vm + i]V^b. (7) 

The "semi-Lagrangian" time derivative G, = d, + Ui • V ac- 
counts for advection by the shearing background. The shear- 
ing component in the Coriolis term cancels out the tidal force. 
We remark that equations (|6]l and (|7]i are not just linearized 
equations, they remain valid even if the amplitude of the per- 
turbations is not small compared to the background values, 
and they are exact as long as a single Fourier mode is con- 
sidered. Under these conditions, it is sensible to study the 
evolution of u(x, t) and b(x, t) for a long time. 

In order to solve equations (|6]l and (|7]i, it is convenient 
to work in Fourier space. The x-dependence of the "semi- 
Lagrangian" time derivative can be removed by employing a 
shearing coordinat e system jx' ,y' , z' , t') = (x, y + q Q.oxt, z, t) 
in which D, = d,- jGoldreich & Lvnden-Beli|[l965l) . A single 
mode with a fixed "shearing" wavenumber k' is thus given by 



u{x, t) - 2Re [uk'{t) exp(/^' • x')] , 
b{x, i) = 2Re {bk'(t) exp(/A;' ■ x')] , 



(8) 
(9) 



where k' ■ x' - k(t) - x - {k'^ + q ^otk'y)x + k'^y + k[z. 

Substituting the ansatz (O and (|9]l into equations (|6]l and (|7]l 
leads a set of equations for the Fourier amplitudes 



d,u - q ClQU^y - icJAb - vk u - 2Q()Z xu - ikP 
d,b + qQ.Qbxy - icj\U - r/k^b. 



(10) 
(11) 



Here, we have replaced d,' by d, and omitted the subscripts 
k' in the Fourier coefficients in order to simplify the notation. 

' Formally speaking we consider two modes, with wavenumbers k and 
-k. However, because the functions under consideration are real, the Fourier 
coefficients satisfy /_t = fl- 



We have also introduced the (tim e-independent) Alfven fr e- 
quency wa = Bo(t) ■ k(t) (see also lBalbus & Hawleylll992al) . 

The pressure term can be eliminated from equation (fTOl i us- 
ing the solenoidal character of the velocity field, which im- 
plies d,ik = ikd, + iqQ.okyX. This leads to 



2/0, 



k^ 



-[(q- l)kyU^ + k^Uy^, 



(12) 



which is independent of u,. Because we are interested in 
the transport of angular momentum along the radial direction, 
this decoupling allows us to solve separately for the x- and 
y-components. 



3. DYNAMICAL EVOLUTION OF MRI-STABLE MODES 

The dynamical evolution of the modes with k = k,z 
is quite simple; they grow exponentially if fe^w^ < 2^Qq 
(IBalbus & HawlevlfT9^ IPessah et alj r2006>. Thus, a Kep- 
lerian disk (with q = 3/2) can exhibit exponential growth but 
a shear profile with ^ < only supports stable oscillations. In 
order to isolate the interesting dynamics that could arise from 
modes that are not associated with the MRI, we thus focus on 
modes with k- - 0. Taking the curl of the momentum equa- 
tion, it is easy to verify that the Coriolis term does not play a 
role in the equation for the vorticity, and hence it has no effect 
on the dynamics of the system (Lithwick 2007). This shows 
explicitly that the standard MRI is absent in our analysis. 

We choose the origin of time so that k^{t) is initially zero. 
In other words, we use k^it) to define our time coordinate. 



T = k,(t)/ky = q Qof , 
so the divergence-less conditions become 

T M + U^, — t Bx + by — 0. 



(13) 



(14) 



Assuming v, // <K LOf^jk^, we can work in the ideal limit and 
neglect viscosity and resistivity. The ;c-components of the 
MHD equations then become 

where all the temporal dependence is contained in the factor 
r(r) = 2r/(r2 + 1), (16) 



and 



oj = cijA/q^o 



(17) 



is the dimensionless Alfven frequency. The linear system ( fTST i 
can be recast into one equation as 



d^b^ + r(r) drb^ + (Jb,, = 0, 



(18) 



where the dependence on the parameters q, Qq, and wa is 
only through the combination {(OaI qQ.o)^ . This second order 
differential equati on for b^ is identical to equation (2.20) in 
Balbus & H awlevI (fl99 2a) when A:. = 0. In this case, the per- 
turbations in the z-coordinate decouple from the perturbations 
in the perpendicular direction, see also their equation (2.19). 

Unfortunately, equation ( fTST i does not have an analytical so- 
lution. However, if we consider the Hmit » 1, it reduces to 
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Figure 1. The Fourier amplitudes for the magnetic field components by (left) and hy (right). The thick solid lines are the numerical solutions obtained by solving 
equation 118) with q = - 3/2 a nd wa_= Hq and using the divergence-less conditions 1141 . The asymptotic solutions for early and late times are both given by the 
analytical approximation j20) and \22\ but with different integration constants. The dotted lines correspond to C- = 1 and 5- = 0, which are used to specify the 
initial conditions alLiir = —20. The dashed lines correspond to C+ = 0.285 and S + = -1.896, which are obtained by matching the analytic and numerical results 
at late times, e.g., u>r = 20. Note that, for the sake of clarity, we truncate the analytical solutions in the right panel at ttir = 2 and —2. 



a spherical Bessel equation, which posseses as solutions 

Ux = SjiidJT) + Cy]{a)T), 

C + Sa)T S - CcDT . 
= 7—5 — cos(coT) + — — sin(wT), (19) 



bjc = -iSjoiair) - iCyoicDT) 
- — cos(wt) sin(w r), 

(l)T COT 



(20) 



where j„(x) and y„(x), with n - 0, 1, are spherical Bessel 
functions of the first and second kind, respectively; and S and 
C are complex constants determined by the initial conditions. 
Using the divergence-less conditions in equation (fl4] i. the y- 
components are simply 

Uy = -T[Sjl((DT) + Cyi(OJT)] 

C + SOJT , ^ S - CCJT . , ^ , ^ 

cos(lot) sm(wT), (21) 

by = IT [S t) + Cyoio) t)] 

iC iS 
— cos(tij t) -I sin(w t). (22) 

Because the pressure in equation ( fT2l i is independent of u,, the 
exact solutions (for all time) are Alfven waves. 



M. = C'cos(a)T) + 5'sin(wT), 
b~ = iC sin(faj t) - iS ' cos(w r). 



(23) 
(24) 



where C and 5" are some other complex constants. We note 
that, although the direction of the dimensionless time r de- 
pends on the sign of the shear parameter q, the combination 
ix)T = oJht is insensitive to it. Hence, given the same initial 
conditions, and I? v are symmetric, while Uy and by are anti- 
symmetric, in the shear parameter q. 

We demonstrate the accuracy of these analytical approxi- 
mations in Figure [1] which shows both the numerical and an- 
alytical solutions for Im[l7v] and Im[/7v] with q - -3/2 and 
wa - as an example. The initial conditions are set at 
(jL>T - a)\t - -20 by choosing - 1 and S - - 0. The nu- 
merical solutions, shown with thick solid lines, result from in- 
tegrating equation (fTSl l with the definition ( fT6] l and using the 
divergence-less conditions ST4i . The dotted lines in the two 
panels are obtained by setting C = C_ and S = 5 _ in the an- 



alytical approximations (l20l i and (l22l i. These solutions are in- 
distinguishable for T < -1. As expected, the approximations 
break down for r ^ 0. This is precisely where the numerical 
solutions change their amplitudes significantly. The analyti- 
cal expressions ( l20t and (l22l i are again in excellent agreement 
with the numerical solutions for t > 1, provided that their am- 
plitudes are given by C = C+ = 0.285 and S - S + = -1.896. 
These constants are found by requiring that both the numer- 
ical and analytical solutions match for w r = ojaI » 1 (in 
practice we set = 20). 

Even though our analytical approach cannot predict the 
change in amplitude close to t « 0, the solutions that we ob- 
tain are a very good approximation to the numerical results as 
long as > 1 . We could in principle obtain the coefficients 
C+ and 5 + by an asymptotic matching techn i que si milar to the 
one employed in lHeinemann & Papaloizoul(l2009h . However, 
the analytical solution near t ^ contains special functions 
that are too complicated to be useful. More importantly, as 
we show below, the most interesting features of the solutions 
are independent of the precise values of these constants. 

4. LATE TIME STRESS AND ENERGY 

Given the solutions (fT9] l - (l24l l for the Fourier amplitudes, 
we obtain the (mean) total stress Tj^-y = {u^Uy - bjj^) and 
(mean) energy density E = {u^ +b^}/2 of the fluctuating 
field fl where the brackets stand for the spatial average, see, 
e.g., IPessah et al.l (120061) . Because these solutions are only 
valid for early /late times, we can approximate the total stress 
and energy density up to first order in 1 / w t as 

T,y ~ — ^ [(|5 P - left cos(2w t) + (S'C + S C*) sin(2w t)1, 

(25) 

£ ~ — \-\i\S\^ -\C\^) sin(2w t)-(S*C + SC*) cos(2w t)1 

CO^T^ J 

+ \(\S\'-\Cf) + \Sf + \C'\\ (26) 

where the asterisk denotes complex conjugation. Using these 
expressions, it is easy to see that the energy balance equation 
d,E - qQ.oTyy is also satisfied up to order I/ojt. 

^ We do not include the stress and energy generated by the time-dependent 
mean field Bo{t) here. This contribution depends on the initial magnetic field. 
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Figure 2. The thick and thin solid lines correspond, respectively, to the total 
energy E{t) and total stress T„(t), calculated using the Fourier amplitudes 
fij-, iiy, h^, and by obtained numerically. The dotted and dashed lines show 
the analytical approximation for the energy, equation )26K using the two sets 
of constants described in the caption of Figure [T] The ,t-components of the 
velocity and magnetic fields are symmetric in q but the y-components are anti- 
symmetric. Therefore, changing the sign of the shear parameter q changes the 
sign of the stress T„ but not of the energy E. 

In Figure |2] we illustrate the numerical solutions for the 
stress Tsvit) and energy E{t), given by the thin and thick solid 
lines, together with the analytical approximation for the en- 
ergy. The latter has been obtained by substituting the two 
pairs of constants, C- - I and 5'_ = 0, and C+ = 0.285 and 
S+ = -1.896, into equation ( |26] |. It is thus clear that the late- 
time stress oscillates around zero with decreasing amplitude, 
while the energy density asymptotes to a non-vanishing, time- 
independent value. The expression for the energy density at 
early/late times in terms of the constants S ± and C± is given 



= lim Eit) 

r— >±oo 



\s±\' + \cj 



(27) 



Therefore, the energy gain via swing amplification, is 
in general a function of the ratio w = coA/q^o and the initial 
conditions. However, it is possible to obtain conclusions that 
are independent of the latter 

The dependence of the energy gain on the initial condi- 
tions for aP' - 1 is shown in Figure [3] The horizontal axis 
describes how the initial energy is distributed between the 
j„ modes and the y„ modes; while the different lines show 
the phase difference in the corresponding initial amplitudes. 
When arg(5'-/C-) = n/2 or 3n/2, the j„ and y„ modes are 
completely out of phase and evolve independently. This re- 
sults in an energy gain which is linear in the initial amplitudes 
(thick solid line). We have found that the dependence of the 
energy gain on the phase difference between the constants de- 
termining the initial conditions is weaker if ll) decreases below 
unity. In this case, all the different curves converge to the thick 
line corresponding to arg(5'_/C_) = n/2. At the same time, 
as 0) decreases below unity, this line gets steeper, providing 
thus a larger energy gain. This justifies referring to the y„ and 
j„ as the "growing" and "decaying" modes, respectively. 

We illustrate the dependence of the energy gain on the shear 
parameter in Figure|4](because the results depend only on u/, 
we only show the positive domain in the horizontal axis). In 
the limit of weak shear, there are only pure Alfven waves and 




14 0.6 

SV(S-+C'_) 

Figure 3. The solid curves show the dependence of the energy gain, E+/E-, 
on the initial conditions for oj = aix/q^k) = 1- The horizontal axis provides 
a measure of the relative amplitude of C- and S - ; while the dilferent fines 
sfiow the results for various phase dift'erences. For strong shear 1 /a/ > 1 , the 
dependence of the phase is weaker and all the dift'erent curves collapse onto 
the straight thick fine, which corresponds to arg(5-/C-) = !t/2, while this 
one gets steeper, thus providing a larger energy gain, see also Figure|4] 

there is no net energy gain. The dashed line shows that the 
energy gain tends to the value E+/E- - lO/w^ as 1/w » 1, 
which provides a good description of the numerical results 
for strong shear. The asymptotic behavior is insensitive to the 
initial conditions as long as the growing mode is excited, i.e. 
C + 0. 

5. DISCUSSION 

5.1. Summary and Connection to Previous Work 

We have employed the shearing-sheet framework to study 
the dynamical evolution of MHD waves in weakly magne- 
tized differentially rotating backgrounds which are stable to 
the MRI. While the fact that these waves can be transiently 
amplified is widely appreciated, our motivation to study them, 
as well as the results that we obtained, concern dynamical as- 
pects that have not received as much attention. This is whether 
these shearing MHD waves can play a significant role in the 
transport of angular momentum in regions of the disk where 
the MRI is inefficient, such as the accretion disk boundary 
layer 

The eq uations that we hav e solved are similar to those pre- 
sented in IBalbus & HawlevI (I1992al) . who provided numer- 
ical solutions showing that transient amplification of MHD 
waves is a general outcome for the modes with wavevectors 
that are not exactly aligned with the rotation z-axifl Ana- 
lytical insight on the non-linear dynamics of these waves has 
been usually hindered by the fact that the governing equa- 
tions cannot be simpl ified beyond a se t of coupled differential 
equations (see, e.g.. Fan & Lou| [T997l: IKim & Ostrike:^ l2000l: 
Brandenburg & Dintrans 2006, and also]johnson 2007, where 
higher order WKB solutions for the linear evolution of the 
shearing waves are provided). In this paper, by isolating the 
modes that are unrelated to the standard MRI, i.e., by setting 
A:. = 0, we have been able to provide analytical solutions that 
are valid for all times, except close to the instant where the 
waves evolve from leading to trailing, i.e., when kjc(t) = 0. 
These s olutions are exact when only one Fourier mode is con- 
sidered (IGoodman & Xul[T994l) . 



Here, we assume S ' = C = 0, and thus avoid the uninteresting contribu- 
tions due to any stable Alfven wave initially present, see equations j23M 24t. 



For an analysis of non-axisymmetric spiral wa ves when only a str ong, 
vertical background magnetic field is considered see llagger et alj (199'2\) . 
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Figure 4. The filled circles represent the value of the energy gain E+ /£_ 
for different values of the shear, parameterized via qClo/aix = 1/lo, obtained 
via numerical integration using the initial conditions C- = 1 and S- = 0, 
i.e., only the growing mode is excited. In the limit of weak shear, there are 
only pure Alfven waves and thus there is no net energy gain. The dashed line 
shows the function 10 /co^ = \0(qSli:,l(jjj^)^ , which is in good agreement with 
the numerical results for strong shear. This asymptotic behavior is indepen- 
dent of the initial conditions as long as C- 4^ 0. 

Despite the fact that we do not predict analytically the 
amplification factor for these waves, the characteristics of 
the solutions that we found allowed us to draw important 
conclusions. We showed that the amplification factor is 
only a function of the (time-independent) dimensionless ratio 
(q'Qo/wa)^, with 



E- \ WA / 



for ^Qo » (jDa , 



(28) 



see equation ( |27l ) and Figure|4] and it is thus insensitive to the 
sign of the shear parameter q. 

An important result of this study is that while the energy of 
these MHD waves can be significantly amplified, their net as- 
sociated stresses oscillate around zero, see equation dZSl l and 
Figure |2] This suggests that these shearing MHD waves are 
unlikely to play an important role in the transport of angular 
momentum in the accretion disk boundary layer region. These 
findings are consistent with the results of global MHD simu- 
lations ofaccretion_disks w ith a rigid inner boundary carrie d 
out in lArmitagd (I2002h and ISteinacker & Papaloizoul (I200a . 
These simulations show that the inner disk regions, where 
dQ./dr > 0, can develop strong toroidal magnetic fields, with 
associated magnetic energies that can easily reach a few tenths 
of the thermal energy, without leading to efficient angular mo- 
mentum transport. 

5.2. Implications 

The importance of understanding the relationship between 
the stress and the radial gradient in angular frequency re- 
sides in that this dependence plays a key role when mod- 
eling the inner structure of an ac cretion disk surrounding a 
weekly magnetized star (see, e.g.. IPopham & Naravanlll995t 
IPopham et al][T996h . 

In the steady state, the constant inward flux of angular mo- 
mentum at any given radius ro is given by 7 = Ml-lnr^HTnp, 
where M stands for the accretion rate, / is the specific an- 
gular momentum, H is the disk height, and r,^ (denoted by 
Txy in our analysis) is the component of the stress respon- 
sible for the flux of azimuthal momentum across the radial 



direction. Thus, the angular momentum flux has two contri- 
butions: 7matter = Ml, which accounts for the flux of angular 
momentum due to mass accretion, and 7sti-e,ss - -^TTr^HTr^, 
which accounts for the flux of angular momentum due to the 
stress acting on the fluid elements constituting the disk. Un- 
der the reasonable assumption that the disk must be Keplerian 
well beyond the boundary layer, i.e., for r » r^, and that the 
stress should vanish in the absence of shear, we must have 
j s Mlifb) > 0, i.e., a slowly rotating star accretes mass and 
angular momentum. 

In order for this picture to be self-consistent, the stress 
must satisfy Tf^{ri, < r < r\,) < 0. In the standard ac- 
cretion disk model this requirement is satisfied by assum- 
ing that the stress is linearly proportional to the local shear 
r,0 ~ -dQ./dr. With this model for the stress, and some 
supplementary assumptions, it is possible to solve for the ra- 
dial dependence of Q(r) and determine the structure of the 
disk (see, e.g., Popham & Naravan 1991). However, this 
assumption, broadly adopted in the framework of enhanced 
turbulent disk viscosity, does not seem to be supported by 
the modern paradigm, in which angular momentum trans- 
port is due to MHD turbulence driven by the MRI. Indeed, 
both local numerical simulatio ns of shearing-boxes with non- 
Keple rian shear profiles (.Pessah et al.l I2008t ISnellman et alJ 
120091) and globa l disk simulations with a ri gid in ner boundary 
(IArmitagd[2002L ISteinacker & Papaloizoull2002l) suggest that 
angular momentum transport is inefficient if dQ/dr > 0. This 
suggests that the detailed structure of accretion disk boundary 
layers resulting from the interaction of an MHD disk with a 
weakly magnetized star could differ appreciably from those 
derived within the standard turbulent shear viscosity, where 
the direction of angular momentum transport is always oppo- 
site to the angular frequency gradient. 

5.3. Final Remarks 

It is worth mentioning explicitly that the shearing-sheet 
framework that we have employed is inherently limited to 
address the conditions expected in the accretion disk bound- 
ary layer For example, the absence of a hard-inner boundary 
could pre vent Kelvin-H elmholtz instabilities from operating, 
see, e.g., IBalsara et al.l (12009). However, these instabilities 
do not seem t o play a predominant role in the gl obal MHD 
simulations of lArmita ge (2002) and Stei nacker & Papaloizoul 
(|2002||) . Moreover, because we have assumed a constant 
background density, our analysis precludes the possibility of 
buoyant modes or convective instabilities. Whether these 
instabilities, and the turbulence they could drive, transport 
angular momentum inward or outward in Keplerian disks 
has been lo ng debated (Rvu & Goodman 1992; Cabot. 19961: 
IStone & Bal bus 1996; Lesur & Ogilvie 2010). To our knowl- 
edge, these convective instabilities have not been studied in 
differentially rotating backgrounds with angular frequencies 
increasing outward; and speculating about their role goes be- 
yond the scope of the present work. 

In spite of the simplifications of our analytical study, the 
explicit solutions that we have found can provide physical in- 
sight and help elucidate transport processes in the inner disk 
regions close to a weakly magnetized accreting st ar. The cur- 
rent availability of powerful parallel codes (e.g. IStone et aLl 
I2008ii) which are already being used to stu dy the hydrody- 
namic s of accretion disk boundary layers (Be lvaev & Rafikovl 
1201 11) holds the promise that a more detailed understanding of 
MHD boundary layers will soon be possible. 
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